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ABSTRACT 

We report on follow-up observations of the gamma-ray burst GRB 060927 
using the robotic ROTSE-IIIa telescope and a suite of larger aperture ground- 
based telescopes. An,., optical afterglow was detected 20 s after the burst, t-he 
earliest rest-frame detection of optical emission from any GRB. Spectroscopy 
performed with the VLT about 13 hours after the trigger shows a continuum break 
at A ~ 8070 A, produced by neutral hydrogen absorption at z ~ 5.6. We also 
detect an absorption line at 8158 A which we interpret as Si II A 1260 at z = 5.467. 
Hence, GRB 060927 is the second most distant GRB with a spectroscopically 
measured redshift. The shape of the red wing of the spectral break can be 
fitted by a damped Lya profile with a column density with logfAGi/cm" 2 ) = 
22.50 ± 0.15. We discuss the implications of this work for the use of GRBs as 
probes of the end of the dark ages and draw three main conclusions: i) GRB 
afterglows originating from z > 6 should be relatively easy to detect from the 
ground, but rapid near-infrared monitoring is necessary to ensure that they are 
found; ii) The presence of large H I column densities in some GRBs host galaxies 
at z > 5 makes the use of GRBs to probe the reionization epoch via spectroscopy 
of the red damping wing challenging; Hi) GRBs appear crucial to locate typical 
star-forming galaxies at z > 5 and therefore the type of galaxies responsible for 
the reionization of the universe. 

Subject headings: gamma rays: bursts (GRB 060927) — cosmology 
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1. INTRODUCTION 

It is well established that most long-duration gamma-ray bursts (GRBs) are caused by 
the death of massive stars (e.g. Galama et al. 1998; Hjorth et al. 2003; Stanek et al. 2003) 
and, due to their brightness, they can be seen throughout the observable universe (Kawai et 
al. 2006). Given these facts, it has long been realized that GRBs could be powerful probes 
of star-formation activity throughout the history of the universe (e.g. Wijers et al. 1998). 
The currently operating Swift satellite (Gehrels et al. 2004) has revolutionized the GRB 
field: it has not only increased the detection rate of rapidly well- localized GRBs by roughly 
an order of magnitude compared to previous missions, but it also detects much fainter and 
more distant bursts (Jakobsson et al. 2006a). In the Swift era, GRBs have indeed become 
cosmological probes of the early Universe, as already extensively predicted in the literature 
(e.g. Wijers et al. 1998; Lamb & Reichart 2000). In particular, Barkana & Loeb (2004) argue 
that GRBs may be optimal probes of the epoch of reionization. 

Here we present observations of GRB 060927 for which we determine a very high redshift 
(z = 5.467, § 2.3). In § 3 we discuss the light curve and the properties of the afterglow in 
the context of the fireball model. In § 4 we end by revisiting the question of how GRBs may 
be used to probe the end of the dark ages. 

Throughout this paper we assume a cosmology with H 0 — 70 km s~ 1 Mpc -1 , il m =0.3 
and fi A = 0.7. In this model, a redshift of 5.47 corresponds to a luminosity distance D\ nm = 
51.9 Gpc and a distance modulus jj, = 48.6 mag. At that distance, 1" on the sky corresponds 
to 6.01 proper kpc and the look-back time is 12.4 Gyr (roughly 92% of the time since the 
Big Bang). 


2. OBSERVATIONS AND DATA ANALYSIS 
2.1. High-Energy Properties 

GRB 060927 and its X-ray afterglow were detected by the Burst Alert Telescope (BAT) 
and the X-ray Telescope (XRT) on-board the Swift spacecraft on 2006 September 27.58860 
UT. At the measured redshift of z = 5.47 (§ 2.3) the BAT 15-150 keV energy band and 
the XRT 0.3-10 keV energy band correspond to 97-970 keV and 1.9-64.7 keV in the source 
rest frame, respectively. Here we describe the burst phenomenology as observed by the Swift 
instruments. Data have been analyzed using the standard analysis software distributed 
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within FTOOLS 1 v. 6.1.1. 

The gamma-ray prompt emission shows a complex light curve with an initial bright 
episode, split into two peaks. The peak flux is reached during the first peak, ~ 0.7 s after 
the trigger (To), and the emission subsequently decays to the background level at ~ T 0 + 9 s. 
A fainter hump is visible at ^ To + 20 s. The observed duration is T 90 = 22.6 ± 0.3 s, which, 
accounting for the (1 + z) time stretching factor, corresponds to 3.54 ± 0.05 s in the source 
rest frame. GRB 060927 is thus a long-duration GRB (Kouveliotou et al. 1993), both in the 
observer’s frame and in the rest frame. 

A simple power-law model is not a good fit to the BAT time-average spectrum of the 
burst (y 2 = 72 for 57 d.o.f.). A significant improvement is achieved by adopting a cutoff 
power law (y 2 = 58 for 56 d.o.f.), which provides a photon index T = 0.9 ± 0.4 and a peak 
energy E p = 72 + g 9 keV. The observed burst fluence is 1.1 x 10 -6 erg cm -2 in the 15- 
150 keV band. A fit with a Band model (Band et al. 1993) does not significantly improve 
the y 2 value, albeit this can be due to the low sensitivity of BAT at high energies. Using 
the best-fit model, we can classify GRB 060927 as an X-ray rich GRB using the definition 
proposed by Lamb, Donaghy, & Graziani (2005). During the initial prompt emission, hard- 
to-soft spectral evolution is present, the first spike being significantly harder than the second. 
The spectrum of the late bump is consistent with that of the second spike. 

The X-ray afterglow was detected by XRT about 70 s after the trigger and it was mon- 
itored for the following 3 days. The X-ray afterglow position, calculated using the updated 
boresight (Burrows et al. 2006), is R.A. (J2000) = 21 h 58 m 12^03, deck (J2000) = 05°21 , 49'.'4, 
with a 90% error radius of 3"7. This position is within O'.'? of the optical counterpart (§ 2.2), 
The X-ray light curve displays an initial shallow decay with a slope cti = 0.71 ± 0.06, steep- 
ening to a 2 = 1.35 ± 0.11 at ~ 4 ks after the trigger. The spectrum is modeled with an 
absorbed power law with photon index T = 1.87 ± 0.17. After correcting for the Galactic 
absorption N u = 5.2 x 10 20 cm' 2 , the 90% confidence limit for the intrinsic absorption is 
A'h < 3.4 x 10 20 cm -2 (observer frame). Indeed, little observed X-ray absorption is ex- 
pected at this large redshift (Grupe et al. 2006; although see GRB 050904: Watson et al. 
2006; Cusumano et al. 2006; Campana et al. 2007b), since most of the important absorption 
edges have been redshifted out of the observed X-ray region. The unabsorbed flux in the 
0.3-10 keV energy band from To -I- 75 s to To + 11.3 ks is (5.7 ± 0.8) x 10 -12 erg cm -2 s -1 
from which we derive a lower limit of 1 0 -6 erg cm -2 to the fluence emitted by the X-ray 
afterglow. 


1 http : / /heasarc . gsf c .nasa . gov . 
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2.2. Ground-Based Imaging Observations 

The optical afterglow was first detected by the 45 cm telescope of the Robotic Optical 
Transient Search Experiment in Australia (ROTSE-IIIa). The first unfiltered image started 
16.5 s after the burst trigger (Schaefer, Yost, & Yuan 2006), which corresponds to 3 s in the 
GRB rest frame. Therefore, this is the earliest rest-frame detection of optical emission from 
any GRB 2 . 

Subsequent follow-up observations were obtained by several ground-based facilities. In 
particular, we used the data from the following instruments: the robotic 2.0 m Faulkes 
Telescope South (FTS) located at the Siding Spring Observatory, Australia, equipped with 
RatCam; the 105 cm Schmidt Telescope at the Kiso Observatory in Japan; the 2.56 m 
Nordic Optical Telescope (NOT), equipped with ALFOSC at La Palma (Canary Islands); 
the 4.2 m William Herschel Telescope (WHT), using the Aux port imager, also at La Palma; 
the 1.54 m Danish Telescope situated at La Silla Observatory in Chile; and the 8.2 m Antu 
and Kueyen Unit Telescopes of the ESO Very Large Telescope (VLT) on Cerro Paranal in 
Chile, equipped with FORS2 and FORSl, respectively. Near-infrared (NIR) observations 
were secured using the ISAAC camera on the ESO VLT Antu. In Table 1, we provide the 
complete observing log. The reduction of the optical and NIR data was performed using 
standard IRAF 3 tools. The photometry of the afterglow was performed using both PSF- 
matched and aperture photometry. 

Using the VLT image from 2006 September 30, we determined the position of the 
afterglow to be R.A. (J2000) = 21 h 58 m 12!02, decf. (J2000) = +05°2T48'.'9, relative to « 30 
isolated, non saturated USNO-B1 stars. The estimated uncertainty is about O'/ 2 in each 
coordinate. 

The detection of the afterglow by ROTSE and the lack of an optical counterpart ob- 
servation by UVOT on board Swift (Oates & Barbier 2006a) revealed within a few minutes 
that the afterglow of GRB 060927 was very red. This was soon confirmed by the FTS obser- 
vations, which showed an /-band afterglow with no corresponding R-band emission within 
the first 30 minutes after the GRB (Guidorzi et al. 2006a). Later, deep WHT observations 
taken 8.9 hr after the burst indicated R — / >2.0. The R-band image from the NOT taken 
6.2 hr after the burst reveals a low signal-to-noise detection of the afterglow with R ~ 24.4. 


2 See http : / /gr ad.40 . as . ut exas . edu/ grbl og . php . 

3 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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After realizing the high redshift nature of the GRB (see below) we obtained deep I-band 
imaging with the VLT leading to a detection of the afterglow ~ 2.6 days after the burst, and 
an upper limit at 15 days. We also obtained VLT/ISAAC NIR imaging leading to tentative 
detections in the J and K bands. No host galaxy was detected down to I > 25.8 (3 a limit). 
In Fig. 1 we show the afterglow of GRB 060927 as detected in the ROTSE and /-band VLT 
exposures. 


2.3. Spectroscopic Observations and Redshift Determination 

Spectroscopic observations were performed with the FORSl instrument at the VLT 
12.5 hr after the trigger. We used the grism 300V covering the spectral range from 3600 
to 8900 A at a resolution of 12 A (using a l'.'O wide slit). The acquisition was performed 
by aligning the slit with a galaxy 3"1 North of the afterglow and using a position angle for 
which the slit covered the afterglow position (PA = 20° West of North). The afterglow was 
not detected in the A-ba.nd acquisition image down to a limiting magnitude R ~ 24.0. 

In the spectrum, a faint continuum is detected at the expected position of the afterglow, 
but only redward of A = 8070 A. This feature is similar to what is seen in the spectra of 
quasars at redshift z > 5 (Cool et al. 2006a). From this clear and sharp discontinuity in 
the spectrum we conclude that GRB 060927 is a high-redshift event, and is not red because 
of dust extinction. We interpret the absorption blueward of the break, where the flux level 
is consistent with zero, as a Lya forest trough at redshift z ~ 5.6. In the sky-subtracted 
two-dimensional spectrum we also clearly detect a single absorption line at A = 8148 A in a 
wavelength region clear of bright sky lines (see Fig. 2). This is most likely Si II A 1260 at a 
redshift of z — 5.467. At this redshift we would also expect O I/Si II at 1303 A and C II at 
1335 A, but these fall on bright sky lines and hence would not be detectable. This value for 
the redshift supersedes our preliminary determination (Fynbo et al. 2006). The difference 
in the redshift derived using the hydrogen and silicon lines is due to the presence of strong 
neutral hydrogen absorption (see § 2.4): the damping wing of the Lya trough effectively 
shifts the cutoff at redder wavelengths, mimicking the effect of a larger redshift. In Fig. 2 we 
show the one-dimensional flux- and wavelength-calibrated spectrum (upper panel) and the 
processed CCD image (lower panel) . In the latter, we can see the sharp continuum break in 
the afterglow spectrum. The two traces above correspond to nearby objects falling on the 
slit. We last note that the pseudo-redshift, for this burst (Pelangeon & Atteia 2006; Atteia 
2003) is z — 2.37 ± 0.75, significantly different (4cr) from the spectroscopic value. 
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2.4. Fitting the Red Damping Wing 

Measuring the shape of the absorption profile of the damping wing provides a determi- 
nation of the column density of the neutral hydrogen along the line of sight to the progenitor 
in the host (Madau & Rees 2000). Using the Fit-Lyman program within MIDAS, we fit- 
ted the wing shape with a Damped Lya absorber system centered at 7861 A, assuming 
z-Dhk = zsm = 5.467. In Fig. 3 we show the fit to the Lya and Si II A 1260 absorption line 
profiles in a ±30,000 km s -1 velocity interval (solid curve). The best- fit column density is 
estimated to have log(A r Hi/cm -2 ) = 22.50 ± 0.15. The Nhi value is firmly set by the core 
of the DLA line, i.e., the part of the spectrum which has zero residual flux, whatever the 
continuum normalisation is. while the error is mostly due to the continuum placement. The 
A'hi value is in the range of those found for other Swift GRBs at z > 2 (Jakobsson et al. 
2006b), although at the high end of the distribution. 


3. AFTERGLOW PROPERTIES 

3.1. Light Curve 

A major complication with the analysis of the optical light curve is the photometric 
calibration. The ROTSE data were taken with no filter. For the other telescopes, the I- 
band data were obtained using three different filter variants: an SDSS-like i filter for the 
FTS, the cutoff Harris I for the WHT, and the “standard” Bessell I for the VLT. For all 
of them, a significant fraction of the light within the effective filter passband is absorbed by 

the Lyct forest. This makes it difficult to obtain a meaningful calibrated magnitude or flux 

from the measured counts. 

The ROTSE telescope spectral response is determined by the CCD characteristics, which 
results in an equivalent broadband filter with a spectral window from 3000 to 10,000 A and 
a peak efficiency around 6000 A. Due to absorption by the Lya forest, we expect little flux 
blueward of 8000 A and therefore our effective response for the optical transient is closer 
to the I band rather than to the R band. To obtain a consistently calibrated afterglow 
flux, we have used the SDSS photometry of neighboring stars (Cool et al. 2006b; Adelman- 
McCarthy et al. 2006) to establish the numerical conversion factors from ROTSE magnitudes 
to photon count rates. Six bright, unsaturated and isolated SDSS stars were selected near 
the GRB for this calibration, spanning a broad range of stellar color. By using the four- 
band griz SDSS data for these stars, we estimated the total flux in the ROTSE passband. 
After convolving these total fluxes with the ROTSE spectral response function a model flux 
for each of the six comparison stars was produced. The RMS variation of these computed 
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values to the ROTSE measurements was 0.09 mag, which, gives a reasonable measure of 
the absolute accuracy of this procedure. We note, however, that since the counterpart was 
close to detection threshold, statistical uncertainties were always dominating the total error 
budget. 

The next step requires some assumptions about the spectral energy distribution of the 
GRB afterglow. We adopted a “typical” power-law form, F v (v) oc , with 3 = 0.75 (see 
also § 3.3). This spectrum was folded first with the Lya absorption in the intergalactic 
medium (IGM) using the model described by Meiksin (2005) and, second, with the ROTSE 
spectral response. In the ROTSE passband, 83% of the GRB flux is absorbed by the IGM. 
Since optical observations must be compared over four different optical passbands, it is 
desirable to pick a reference wavelength A 0 common to all of them in order to minimize the 
effects of the uncertainty in /?. A sensible choice for Ao requires it to be higher than the Lya 
cutoff frequency (corresponding to 8000 A) and lower than the long-wavelength cutoff of any 
of the various I-band filters (i.e., 8650 A for FTS i). From the above constraints, we chose 
arbitrarily Ao = 8190 A. Rather than quoting magnitudes based on large extrapolations to 
actual filter bandwidths, we computed the spectral flux density at the value of A 0 indicated 
above. We also checked the effect of varying the afterglow spectral index j3 by ±0.25. Since 
the observed spectral region is narrowly confined by the Lya cutoff and by the edge of 
the filters, the corresponding fluxes are not much changed (by ~ 1%). The inferred flux 
numerical values are listed in Table 1, and have been corrected for the Galactic extinction 
E(B — V) = 0.06 mag (Schlegel, Finkbeiner, & Davis 1998) as well. 

A similar, but simpler procedure was carried out for the data obtained from the other 
telescopes. These were in fact equipped with narrow-band filters, so the fraction of light 
inside the filter lost due to the intergalactic medium absorption was computed. Taking into 
account the appropriate filter sensitivity functions and the detector quant um efficiencies, 
this fraction is 68%, 47%, and 56% for the FTS, WHT, and VLT, respectively. 

We also note that the afterglow was marginally detected in the R band by the NOT and 
by the Danish 1.54 m telescopes. This detection is very likely due to the extension of the 
R-band filter response redward of 8000 A. Using the tabulated filter efficiency, we estimate 
that ~ 6% of the light is transmitted (3 mag suppression). Considering the observed color 
R — / = 2.3 ±0.3, and correcting for the Lya suppression in both filters, we infer an intrinsic 
R—I = 0.2±0.3, which is consistent with normal afterglow colors given the large uncertainty. 
We thus confirm that very small light is transmitted blueward of the Lya cutoff. 

In Fig. 4 we present the light curve of the GRB 060927 afterglow, showing /-band 
corrected fluxes and X-ray data points. The X-ray afterglow seems to be rather typical 
(Zhang et al. 2006; O’Brien et al. 2006; Nousek et al. 2006). At, early times ( t < 3000s) 
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it is characterized by a shallow decay, with temporal index aq tv 0.7. Comparison with 
the extrapolation of the BAT flux also implies a steep decline before the beginning of the 
observation, showing that GRB 060927 had all the three “canonical” phases. The shallow 
decline steepened to a 2 ~ 1.3 at ~ 4 ks after the trigger. 

The situation is different in the optical. After ~ 500 s the data are fairly well described 
by a single power law with decay index a = 1.17 ± 0.03 (y 2 = 3.6 with 6 d.o.f.). The early 
ROTSE data, however, lie significantly below the extrapolation of the late light curve; a 
single power-law fit to the full I-band data set provides an unacceptable x 2 = 86 for 11 
d.o.f.. A flat light curve between 300 and 1000 s is also visible in the Kiso A-band data. The 
last ROTSE points agree well with the FTS data, so that the difference cannot be due to an 
intercalibration problem. 

One possibility is that a break was present also in the optical, at ~ 500 s after the trigger. 
Alternatively, the peak at ~ 300 s after the GRB may suggest a rebrightening. The flux after 
the peak follows a single power law, indicating that no flare was present. Rather, the flux was 
steadily larger than the extrapolation of the behavior before the brightening. This step-like 
behavior is consistent with being due to an episode of energy injection. Energy injection 
is also often invoked to explain the shallow X-ray decay (e.g. Zhang et al. 2006), albeit in 
GRB 060927 the shallow phase ended significantly later than the brightening episode. A 
different possibility is that the early emission was due to another component, like a reverse 
shock or central engine activity (note that the first ROTSE measurements are simultaneous 
to the last peak of the prompt emission). 

It is remarkable that the observed I-band flux is declining by the very early observation 
(3 s in the GRB rest frame). This is difficult to explain in the context of the standard 
afterglow model (e.g. Sari, Piran, & Narayan 1998). The flux is indeed expected to rise 
before the so-called deceleration time Gee- If the early emission was due to the forward 
shock, our data constrain Gee <3s (GRB frame). This translates into a lower limit on the 
fireball initial Lorentz factor To > 1000. Furthermore, a break in the light curve is expected 
when the injection frequency n in j crosses the observed band. The lack of such a break requires 
that £/i n j was already below the optical range at the beginning of our observations. Very small 
equipartition parameters are required for this to happen. These difficulties may also support 
the idea that the early emission was not due to the forward shock. 

The optical and X-ray decay slopes at late times are marginally consistent. This is in 
agreement with the broad-band spectral analysis carried out at t ~ 5700 s (§ 3.3), which 
shows that the overall SED is well described by a single power law. The hard spectral index 
(3 v 0.7 (Table 2) suggests the cooling frequency was above the X-ray region. However, as 
discussed before, at early times the X-ray emission was flatter, and the optical light curve 
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showed no steepening at the X-ray break time. This behavior has already been noted in a 
number of Swift bursts (e.g. Panaitescu et, al. 2006). 

3.2. Burst Energetics 

The time-averaged spectral energy distribution of the prompt emission is well fitted by 
a cutoff power law with a photon index T = 0.9 ±0.4 and a peak energy E p = 72+g 9 keV. In 
the burst rest frame, the intrinsic peak energy was E' p — 46G~|o° keV. The burst fluence over 
the T 90 interval in the 15-150 keV band was 1.11°; 7 x 10 -6 erg cm -2 . This corresponds to 
an isotropic-equivalent energy E- mo = 7.7Slj) x 10 52 erg in the rest-frame 1-10 000 keV band. 
The values of E' p and E iso make this GRB consistent with the Amati relation (Amati et al. 
2002; Amati 2006). 

We can provide only a lower limit to any light curve break, which corresponds to t b > 
2.6 days (the epoch of our last detection). The lower limits for the jet opening angle and for 
the beaming-corrected energy are •& > 5° and E iet > 2.6 x 10 5 ° erg, respectively. The E' p vs 
£jet correlation (Ghirlanda et al. 2004) predicts for this burst A — 10lf days (Campana et 
al. 2007b), thus GRB 060927 is consistent with this relation, although it does not constrain 
it (see also Ghirlanda et al. 2007). 


3.3. The Afterglow Spectral Energy Distribution 

We have constructed the spectral energy distribution (SED) of the afterglow covering 
from NIR to X-ray wavelengths at 5709 s after the Swift BAT trigger. Both a power law 
(PL) and a broken power law (BPL) have been fitted including the effect of dust absorption, 
adopting the extinction curves of the Milky Way (MW), the Large Magellanic Cloud (LMC) 
and the Small Magellanic Cloud (SMC) as parametrized by Pei (1992), and of photoelectric 
absorption in the X-ray band (assuming Solar metallicity). All fits include a fixed Galactic 
absorption component, with Ah, G al = 5.2 x 10 2 ° cm~ 2 (Dickey & Lockman 1990) and E(B — 
V ) = 0.062 mag. Correction for Lycc blanketing was applied to the I and R data. When 
excluding the R band data, which are severely suppressed, the parameters did not change 
within the errors. All the results are listed in Table 2, and the best fit is shown in Fig. 5. 
Also considering the broad-band SED, we do not measure any intrinsic X-ray absorption, 
consistent with the XRT analysis alone. A broken power law with photon indices i'j and T 2 
provides no improvement upon the single power law fits, and the break energies are difficult 
to constrain. We also test for the presence of a cooling break in the observed energy range 
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by fixing Tj — f ’2 — 0.5, as expected from the standard synchrotron theory, but found that 
no such break is required. The optical absorption is found to be low, in agreement with 
many previous studies of GRB host galaxy extinction (e.g., Kann et al. 2006; Starling et 
al. 2007). However, it is not possible in this case to differentiate between the three tested 
extinction curves. The 2175 A bump characteristic of the MW extinction curve would occur 
at ~ 1.4 pm when redshifted to z ~ 5.47, falling in between the J and K bands where 
we have no observations. From our data we can, however, obtain an accurate measurement 
of the power-law slope at this epoch, which is consistent within the errors with the value 
derived from fitting the X-ray data alone (§ 2.1). The similar decay indices in the optical 
and X-ray ranges also suggests that both components lie on the same spectral segment. 


4. GAMMA-RAY BURSTS AS PROBES OF THE DARK AGES 

Spectroscopy of GRB afterglows has been argued to be a powerful tool to probe the 
red Lya damping wing resulting from the neutral intergalactic medium before reionization 
(Barkana & Loeb 2004). In fact, they have simple power-law continuum spectra and a 
relatively small ionizing effect on their environments (compared to QSOs). This shape al- 
lows one to directly infer the neutral hydrogen column and hence the state of reionization 
(Madau & Rees 2000). However, as pointed out by Totani et al. (2006), this method may 
be compromised by the strong hydrogen absorption due to the dense environments that 
long GRBs usually seem to be situated in (Jakobsson et al. 2006b). Our observation that 
such a high absorbing column is also present in a z > 5 GRB confirms that it is likely 
difficult to use afterglow spectra as a diagnostic of reionization due to the Lya damping 
wing. We also note, however, that GRB-DLA column densities are usually lower than for 
GRB 060927: in the sample by Jakobsson et al. (2006b) only two out of 18 events 10%) 
have Nni > 10 22 - 5 cm -2 . 

The outlook for the use of GRBs as star-formation tracers is much better. For example, 
it has also been suggested by Wijers et al. (1998) that GRBs could be ideal tracers of the 
average star formation density due to their brightness, independently of dust extinction and 
the fact that they originate from a single (or double) stellar progenitor and so do not require 
a detectable host galaxy (see also Conselice et al. 2005). Evidence has been mounting that 
GRBs may only trace star formation in relatively metal-poor environments (e.g. Le Floc’h 
et al. 2003; Fynbo et al. 2003; Tanvir et al. 2004; Sollerman et al. 2005; Fruchter et al. 
2006; Stanek et al. 2006; Modjaz et al. 2007), In the popular collapsar model, GRBs are 
expected to be formed only by stars with metallicity below about 0.2-0. 3 -Zq (Heger et al. 
2003; MacFadyen & Woosley 1999; Hirschi, Meynet, & Maeder 2005; Woosley & Heger 2006; 
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Langer & Norman 2006, but see also Fryer & Heger 2005). However, before z « 5 the 
majority of the star-formation should occur in environments with Z < 0.2Z o (Nuza et al. 
2007), and therefore the redshift distribution of z > 5 GRBs should directly map the star 
formation history of the universe back to the first era of population-II objects, and possibly 
the first population-III stars (e.g. Bromm & Loeb 2006; Belczynski et al. 2006). 

Furthermore, GRBs afterglows pinpoint their host galaxies, provide redshifts for them, 
and in some cases measures of their metallicity, gas density and dynamical state 4 . This 
provides a route to identify and study the types of galaxies that are responsible for the 
bulk of the star formation, and hence for the bulk of the (re)-ionization of the universe. 
This is very important since other methods to locate high redshift galaxies with current 
instrumentation only locate the brightest systems that are too rare to account for the bulk 
of the ionization (Labbe et al. 2006), possibly with the exception of gravitational lensing 
(Stark et al. 2007). 

Swift has been operating for more than two years and we now have a good sample of 
GRBs from which to constrain their fraction at very high redshift. We followed Jakobsson 
et al. (2006a) and selected only Swift GRBs with low Galactic extinction and promptly 
determined XRT positions. In this way we created an unbiased sample discarding bursts 
with adverse observing conditions. By December 2006, there had been a total of 83 of these 
bursts: 46 with spectroscopic redshift measurements (with a mean (z) = 2.31), 18 with 
no optical/NIR afterglow nor host galaxy detection, and 19 with optical/NIR afterglow or 
host galaxy detection but unknown spectroscopic redshift. In Fig. 6 we show the redshift 
distribution for the 83 bursts. The arrows indicate upper limits based on the information 
available: afterglow colors, host galaxy and/or Lya break constraints. The shaded bar 
indicates bursts for which we have no redshift information at all. 

Several bursts in the sample with very red afterglows have shown properties expected for 
high-redshift bursts: GRB 050502B, GRB 050716, GRB 060108, GRB 060602A, GRB 060719, 
GRB 060814, GRB060923A and GRB060923C (Genko et al. 2005; Tanvir et al. 2005; Levan, 
Tanvir, & Fuhrman 2006a; Jensen et al. 2006; Malesani et al. 2006; Levan et al. 2006b; Fox, 
Rau, & Ofek 2006; D’Avanzo et al. 2006, respectively). However, for GRB 060108 (Oates 
et al. 2006b), GRB 060814 (Campana et al. 2007c), and GRB060923A (Levan, Tanvir, & 
Gorosabel 2006c), the detection of the host galaxy in the R band indicates z < 5.5 and hence 


4 Although the existence of the GRB might select enviroments with special properties, recent observations 
(Chen et al. 2006; Watson et al. 2007; Vreeswijk et al. 2007) have shown that the material being probed 
by afterglow spectroscopy is not in the immediate surroundings of the GRB site, implying that it may be 
representative of the host galaxy properties. 
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the red afterglow color is due to dust extinction rather than to high redshift 5 . 

We now try to constrain the fraction / of Swift GRBs at z > 6. From the initial 
sample we removed bursts with measured or constrained redshift, leading to a sub-sample 
of 20 candidates: 18 bursts with no information at all and 2 bursts with very loose redshift 
constraints. We rejected from this set those bursts with duration shorter than 14 seconds 6 7 
and those with (significant) excess column density in the X-ray spectra. As pointed out by 
Grupe et al. (2006), high-redshift GRBs tend to show little or no X-ray absorption, since the 
affected spectral region is redshifted out of the observing window (see however the striking 
counterexample of GRB 050904 at z = 6.295: Watson et al. 2006; Cusumano et al. 2006; 
Campana et al. 2007b). Only six bursts were considered as candidates for being at z > 6. 
Together with GRB 050904, they set an upper limit of / < 7/83 ~ 8%. Do we have any 
evidence that these bursts are at high redshift? We retrieved from the literature optical 
magnitude upper limits 1 and X-ray fluxes 8 (Butler & Kocevski 2007) in order to estimate an 
upper limit to the optical-to-X-ray spectral index pox- High-redshift bursts necessarily have 
Pox < 0.5 (Jakobsson et al. 2004). Unfortunately, the existing observational limits are not 
deep enough to constrain p Q x < 0.5 for any of the candidates. By considering GRB 050904, 
we can thus set a lower limit />1/83~1%, even if we caution, of course, that this limit 
is based on one single event. It has been argued (Butler & Kocevski 2007) that afterglow 
spectra may be intrinsically curved (due to the presence of spectral breaks) and this could 
spuriously mimic the presence of excess column density. If we remove the requirement of 
zero excess absorption, the upper limit on / becomes / < 16/83 ~ 19%, which is of course 
less constraining and very conservative. 

From what is shown above, it is possible that Swift had detected some GRBs with z > 6 
in addition to GRB 050904. How do we improve our follow-up strategies in order to identify 
these bursts and secure spectroscopic redshifts? It is clear from the discussion above that the 
afterglows of such events will be very red due to Lyman-forest blanketing, but that not all 
red afterglows are from high redshift bursts. Fortunately, the dust-extinguished bursts will 
likely be behind large column densities of soft X-ray absorbing material and hence can be 
filtered away (Grupe et al. 2006). The main issue is therefore to detect the NIR afterglows. 

The two most distant GRBs known, GRB 050904 (e.g. Kawai et al. 2006; Haislip et 


5 The afterglow of GRB 060108 was also detected in the B band, implying a low redshift 2 < 3.2 (Oates 
et al. 2006b). 

6 A burst with Tgo — 14 s at z = 6 has a proper time duration of rs 2 s; see also Campana et al. (2007c). 

7 http : / /grad40 . as . utexas . edu/grblog . php . 

8 http : / / astro .berkeley . edu/~nat/swif t . 
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al. 2006; Gendre et al. 2006; Totani et al. 2006; Tagliaferri et al. 2005) and GRB 060927 
(this work), were both detected with relatively small aperture robotic telescopes (TAROT 
and ROTSE-IIIa, respectively). This suggests that we may expect bright afterglows from a 
substantial fraction of GRBs with z > 6, but these will not be detectable with most of the 
currently operating robotic instruments because of the lack of light in the optical window. 
NIR robotic projects (e.g. REM, GROND on the ESO 2.2m, the 1.3m CTIO telescope, 
PAIRITEL, BOOTES-IR) are expected to provide accurate afterglow positions in just a 
few minutes, and, hence, to increase the sample of spectroscopically measured high-redshift 
GRBs. However, three of these are placed within a few hundred kilometers of each other 
in Chile. In order to make the search of events that are only visible in the NIR as efficient 
as in the optical, a more wide-spread net of NIR (semi)-robotic telescopes is required. A 
possible solution, easy to implement on a short time scale, is to equip many of the 1.5-3 m 
class telescopes, that have had less observing demands since the advent of 8-10 m telescopes, 
with red-sensitive deep-depletion CCDs and z-band filters (see e.g. Guidorzi et al. 2006b). 
In this way the probability of detecting bursts out to z ss 7.5 would dramatically increase. 

5. SUMMARY AND CONCLUSIONS 

We have presented X-ray and optical observations of the afterglow of GRB 060927, 
performed with the Swift spacecraft and several ground-based telescopes. The spectroscopic 
observations made with the VLT showed a flux dropout starting at A = 8070 A. We interpret 
the absorption blueward of the break as due to the Lya forest at redshift z ~ 5.6. We also 
detect a Si II A 1260 absorption line at A = 8148 A, corresponding to z = 5.467. 

The modeling of the red damping wing of the Lya absorption found in the optical 
afterglow spectrum leads to a column density with log (Am /cm: 2 ) = 22.50 ± 0.15. This 
is similar to the values observed in GRBs at lower redshift. Such large values make the 
use of the GRBs to probe the reionization epoch via spectroscopy of the red damping wing 
challenging. 

The light curve in the gamma-ray band shows an initial bright episode, split in two 
peaks, with a total duration T 90 = 22.6 — 0.3 s. The X-ray light curve is characterized by a 
shallow decline phase lasting ~ 3000 s, with the temporal index being a 2 0.7. The shallow 
decline steepened to 02 ~ 1.3 at ~ 4 ks after the trigger. The optical data are fairly well 
described by a single power law with decay index a — 1.17 ± 0.03 after ~ 500 s, while an 
extra component is apparent at earlier times. 

We selected a sample of Swift GRBs with low Galactic extinction and promptly de- 
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termined XRT positions. From this unbiased sample, we constrain the fraction / of GRBs 
at 2 > 6. Including GRB 050904, we set a conservative upper limit / < 19%. We also 
point out that GRB 060927 was detected by a 45 cm telescope (ROTSE). This shows that 
high-redshift GRBs can be very bright, and are in principle easy to detect, provided that red 
and near-infrared detectors are employed. .We suggest to exploit the potential of the 2 -band 
filter, capable to detect bursts up to 2 rj 7.5. 
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Table 1. 


Ground-based observations of GRB 060927. Upper limits are at the 3cr level. 


Time (UT) a 

At (s) 

Exposure (s) 

Magnitude 1 * 

F»{ Ao) c 

Filter 

Telescope 

Sep. 27.588794 

16.8 

5 

14.39±0.20 

6800 ± 1300 

none 

0.45-m ROTSE-IIIa 

Sep. 27.588956 

30.8 

5 

14.67±0.22 

5200 ±1100 

none 

0.45-m ROTSE-IIIa 

Sep. 27.589216 

53.3 

15 

15.56±0.19 

2290 ± 410 

none 

0.45-m ROTSE-IIIa 

Sep. 27.589804 

104.1 

25 

17.16±0.41 

520 ± 200 

none 

0.45-m ROTSE-IIIa 

Sep. 27.590624 

174.9 

200 

17.86±0.30 

275 ± 78 

none 

0.45-m ROTSE-IIIa 

Sep. 27.594007 

467.2 

600 

17.07±0.21 

570 ± 110 

none 

0.45-m ROTSE-IIIa 

Sep. 27.602126 

1168.7 

600 

17.99±0.40 

243 ± 92 

none 

0.45-m ROTSE-IIIa 

Sep. 27.654894 

5727.9 

600 

> 17.04 

< 590 

none 

0.45-m ROTSE-IIIa 

Sep. 27.593845 

453.2 

930 

> 20.1 


SDSS R 

2.0-m FTS 

Sep. 27.597834 

797.9 

60 

> 18.4 

< 570 

SDSS i 

2.0-m FTS 

Sep. 27.602725 

1221.3 

120 

19.02±0.17 

320 ± 50 

SDSS i 

2.0-m FTS 

Sep. 27.609716 

1824.4 

180 

19.64±0,22 

180 ± 37 

SDSS i 

2.0-m FTS 

Sep. 27.616009 

2368.2 

120 

19.43±0.27 

220 ± 55 

SDSS i 

2.0-m FTS 

Sep. 27.623009 

2973.0 

180 

20.57i0.47 

77 ± 33 

SDSS i 

2.0-m FTS 

Sep. 27.602650 

1214.0 

90 

19.5±0.3 


R 

1.05-m Kiso 

Sep. 27.627431 

3355.0 

90 

19.4±0.1 


R 

1.05-m Kiso 

Sep. 27.629468 

3531.0 

90 

19.8±0.2 


R 

1.05-m Kiso 

Sep. 27.844120 

2.21 x 10 4 

3x300 

24.4±0.3 


R 

2.5-m NOT 

Sep. 27.957326 

3.18 x 10 4 

4x100 

> 25.5 


B 

4.2-m WHT 

Sep. 27.964109 

3.24 x 10 4 

4x100 

> 24.6 


V 

4.2-m WHT 

Sep. 27.973495 

3.32 x 10 4 

4x100 

> 24.2 


R 

4.2-m WHT 

Sep. 27.978808 

3.37 x 10 4 

4x100 

22.11±0.16 

7.2 ± 1.1 

Harris I 

4.2-m WHT 

Sep. 28.006500 

3.61 x 10 4 

13x900 

24.5±0.26 


Bess ell R 

1.54-m Danish 

Sep. 28.064850 

4.11 x 10 4 

2x60 

>24.0 


Bessell R 

8.2-m VLT 

Sep. 28.108681 

4.49 x 10 4 

3x1800 



G300V 

8.2-m VLT 

Sep. 30.004618 

2.09 X 10 s 

10x30 

21.35db0.29 


K 

8.2-m VLT 

Sep. 30.039178 

2.12 X 10 s 

30x60 

23.15±0.45 


J 

8.2-m VLT 

Sep. 30.159225 

2.22 x 10 5 

21x300 

24.85±0.17 

0.72 ±0.11 

Bessell I 

8.2-m VLT 

Oct. 13.067980 

1.34 x 10 6 

25x300 

> 25.8 

< 0.3 

Bessell I 

8.2-m VLT 


a The reported times correspond to the beginning of the observation. 
b Not corrected for Galactic extinction. 

c Flux density at Ao = 8190 A, measured in /ijy, corrected for Galactic extinction E(B — V) = 0.06 mag. 




Table 2. Model fits to the SED at epoch 5709 s since trigger. All errors are at the 90% 

confidence level. 


Model 8, 


Ext. curve 1 * 

N H C 

E(B - V ) 

Tj 

E h (keV) 

r 2 

X 2 /dof 

PL 


MW 

< 1.2 

U* A ^_0.02 

1.70±0.06 



9.5/20 

PL 


MW 

0 (frozen) 

0-137^5 

1.70 ±0.06 



9.5/21 

PL 


LMC 

< 1.1 

0.087°;°? 

1 i^Q+0.04 



10-7/20 

PL 


LMC 

0 (frozen) 

0 r)R+ 004 
UlUO -0.03 

1 ’ b ^-0,06 



10.7/21 

PL 


SMC 

< 1.0 

0 07+ 0 - 01 
— 0.02 

1,61 7m4 



11.5/20 

PL 


SMC 

0 (frozen) 

0.07±0.03 

1,66 7q;o5 



11.5/21 

BPL 


MW 

< 3.6 

ft 1 1+O.O6 
ull -0.07 

i-«7S2i 

<106 

1.837q' 17 

6.6/18 

bpl (ri = r 2 

1 

© 

bi 

MW 

<4.7 
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Fig. 1. — Top: unfiltered image from ROTSE-IIIa obtained 19 s after the trigger. Bottom: 
VLT /-band images 2.6 days (left) and 15.5 days (right) after the trigger. The circles indicate 
the position of the optical afterglow. North is up, East is left. 
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Fig. 2. — Afterglow spectrum of GRB 060927 taken 12.5 hr after the trigger. The one- 
dimensional spectrum is shown in the upper panel with the identified absorption Lya and 
Si II lines indicated. The CCD image of the spectrum in the corresponding wavelength range 
is shown in the lower panel (bottom trace). 
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Fig. 3. — Afterglow spectrum with the model fit to the Lya and Si II A 1260 absorption line 
profiles (solid curve). The corresponding column density has log(VHi/cm -2 ) = 22.50 ±0.15. 
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Fig. 4. — Optical (8190 A, filled symbols) and X-ray (1.4 keV, crosses) light curves of 
GRB 060927. The early BAT measurements were extrapolated to the XRT range using the 
best-fit spectral model. Optical data were corrected for the Galactic extinction. The dotted 
line shows the best-fit power-law model to the /-band light curve for t > 500 s, while the 
dashed line shows the best-fit broken power-law model to the X-ray light curve. 
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Fig. 5. — Example of fit to the NIR/X-ray SED. The model (thick solid line) shows a single 
power law with Milky- Way extinction. The lower panel shows the residuals of the fit. 
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Fig. 6. — The redshift distribution of Swift GRBs up to 2006 December. The arrows indicate 
upper limits based on the afterglow colors and/or presence of Lya breaks in the spectra 
and/or detection of the host galaxy. The shaded bar indicates the 18 bursts for which we 
have no constraints on the redshift. 


